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ABSTRACT 




A method of analyzing flow through a turbomachine is summarized which is 
suitable for computer programming* The method, which has been reported in NASA 
publications, is based on an equation for the velocity gradient along an arbi- 
trary quasi -orthogonal rather than the normal to the streamline as used in 
previous methods. A quasi -orthogonal is defined to be any curve that intersects 
every streamline between the flow boundaries exactly once, as does an orthogonal 
to the streamlines. With this method a streamline analysis can be made on any 
given stream surface. A quasi-three -dimensional solution can be obtained by 
using the method for a hub-to-shroud analysis, followed by blade -to -blade 
analyses at hub, mean, and shroud. As an example, the method was applied to a 
radial inflow turbine with splitter blades. The complete quasi-three-dimensional 


solution was obtained in a single computer run. 
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INTRODUCTION 


Quasi-three -dimensional methods have been developed and used for analyzing 
flow through mixed-flow turbomachines [l to 4]. The first step in these methods 
is to obtain a two-dimensional solution on an assumed mean stream surface between 
the blades. This two-dimensional solution is based on an equation for the 
velocity gradient along the normal to the projection of the streamlines on a 
plane containing the axis of rotation. This plane is called the meridional plane, 
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and the projections of the streamlines are called meridional streamlines. 

The streamlines and their normals are used to establish a grid for a meridional- 
plane solution. In cases where the distance between hub and shroud is great 
and there is a large change in flow direction within the rotor, the normals 
vary considerably in length and direction during the course of the calculations. 
Therefore, it becomes difficult to obtain a direct solution on the computer with- 
out resorting to intermediate graphical steps. 

The use of normals is not essential to the method, and it appeared pos- 
sible to obtain a direct solution by the use of a set of arbitrary curves from 
hub to shroud instead of streamline normals. These arbitrary curves will be 
hereinafter termed quasi -orthogonals. The quasi -orthogonals are not necessarily 
orthogonal to each streamline but merely intersect every streamline once across 
the width of the passage. The quasi -orthogonals remain fixed regardless of any 
change of streamlines. By using this technique, it appeared possible to develop 
a computer program that would calculate a streamline solution in the meridional 
plane without any intermediate graphical procedures even for turbomachines with 
wide passages and a change in direction from: radial to axial within the rotor 
blade . 

From the meridional solution it is possible to obtain blade surface veloc- 
ities by several methods. One method is by means of a simple approximate formula 
based on the assumption of linear velocity variation between blades and absolute 
irrotational flow £5] . This gives good results within the rotor away from the 
inlet or outlet. Another method is the use of potential flow equations which 
may be solved by relaxation methods [ 6 ], which gives good solutions, but in- 
volve a rather tedious procedure, even with the use of a computer. It was realized 
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that the basic idea used to obtain a meridional solution using quasi -orthogonals 
could be applied to obtain a blade -to -blade solution. In this case, the quasi- 
orthogonals run from blade -to -blade on a stream surface determined by the 
meridional solution. By extending the solution upstream and downstream, a 
good solution throughout the rotor was obtained with a reasonable effort. 

This paper summarizes a quasi-thr.ee-diraensio.nhl analysis method based on 
using quasi -ortho go nals to obtain both the meridional and the blade -to -blade 
solution. Complete details are given in [7 and 8]. The numerical procedure 
required for obtaining solutions using a digital computer are outlined. As a 
numerical example of the application of the analysis method, a high specific 
speed radial-inlet mixed-flow gas turbine with splitter blades is analyzed.. 

NOMENCLATURE 

a parameter, Eq. (ll) 

b parameter, Eq. . (ll) 

c parameter, Eq. (ll) 

Cp specific heat at constant pressure, ft-lb/slug-°R 

g acceleration due to gravity, ft/sec^ 

h static enthalpy, ft -lb/slug 

m distance along meridional streamline, ft 

N number of blades 

An distance between two stream surfaces 

p absolute static pressure, lb/ft 2 

Ap loss in relative total pressure between inlet and any point, lb/ft^ 

q distance along an arbitrary three dimensional curve, ft 

qi,qo intersection of quasi -orthogonal with boundary of flow passage 
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R gas constant, ft-lb/slug-°R 

r radius from axis of rotation, ft 

r-^ radius where mean stream surface is assumed to deviate from mean blade 
shape in radial-flow turbine 
r c radius of curvature of meridional streamline, ft 
T temperature, °R 

t time, sec 

t q blade thickness in the circumferential direction, ft 
V absolute fluid velocity, ft/sec 

W relative fluid velocity, ft/sec 

w mass flow between two stream surfaces, slug/sec 
z axial coordinate, ft 

a angle between meridional streamline and z axis, radians 
p angle between relative velocity vector and meridional plane, radians 
X ratio of specific heats 

Q relative angular coordinate about the axis of rotation, radians 

A prerotation r^V^-^, sq ft /sec 

p mass density, slug/cu ft 

0 ) rotational speed, radians/sec 

Subscripts; 

i inlet 

m component is direction of meridional streamline, see Fig. 2 
n component normal to quasi -orthogonal and An 

r radial component 

z axial component 


d tangential component 
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Superscripts: 

' absolute stagnation conditions 
relative stagnation conditions 

BASIC ASSUMPTIONS AND PROCEDURE 

The basic assumptions are that there is steady relative flow, and that 
the fluid is nonviscous and isentropic. To this a correction for losses is 
made by assuming a loss in relative total pressure varying fron zero at the 
inlet to a maximum at the outlet of the blade passage. The problem of obtain- 
ing an approximate solution to a three-dimensional fluid flow problem is con- 
siderably simplified if certain assumptions are made so that the problem can 
be reduced to the solution of a series of two-dimensional problems on hub-to- 
shroud and blade -to -blade stream surfaces, as illustrated in Fig. 1. The method 
used herein for obtaining a streamline solution on any stream surface, whether 
hub-to-shroud or blade -to -blade, is essentially the same and will be hereinafter 
called the QO Method. 

The initial assumption used here is that the mean surface between the blade 
is a stream surface. With this assumption, a two-dimensional meridional plane 
solution is obtained by the QO Method. This solution defines a blade -to -blade 
surface for each meridional streamline, and then solutions are obtained on these 
surfaces by the QO Method, giving a velocity distribution throughout the passage. 

EQUATIONS FOR STREAM SURFACE ANALYSIS 
If a stream surface (hub-to-shroud, blade -to -blade, or other) is known or 
assumed, together with the streamline spacing normal to the stream surface, the 
problem has been reduced to a two dimensional problem. Presented here is a 
method of solving this two dimensional problem for fairly general cases. The 
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key to the method is an equation for the directional derivative of the relative 
velocity W, which can he derived from Euler's equation for a nonviscous fluid 
in a rotating coordinate system: 


dW (W_ + W ) 2 . -s 

r _ 9 r x 1 dp 

dt r p ot 


d(rWg + oar 2 ) 


dt 

. 

dt 

Here the usual convention is used that 


1 

p M 


1 

p Sz 


df(r ,9,z,t) = df^r(t),0(t),z(t),t*) 
dt dt 


( 1 ) 


where r(t) ,0(t) , z(t) are the parametric equations for the motion of a mass 

particle as a function of time. The coordinate system is shown in Fig. 2. 

If we multiply Eq. (l) hy Wr, W q, and W , respectively, then add and 

z 

combine terms y we get the energy equation 


1 = ^ dr _ 1 d£ 

2 ^.2 dt p dt 

With the assumption of steady isentropic flow, so that i 

p dt dt' 

Eq. (2) can be inte'grated along a streamline from the inlet to obtain 

W 2 - W? = oo 2 (r 2 - rf) - 2(h - h ± ) 

Now, use the fact that 


V 2 = W 2 + ZVgxr - oo 2 r 2 


to get 


V? 

h i " h i - f = 


2 2 2 
Wi + 2coX - oo rq 


( 2 ) 

(3) 

(4) 


2 


( 5 ) 
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Substituting this in Eq. (3) gives 


h = h^ 


cuA + 


(cur)' 




( 6 ) 


Let q denote the distance along an arbitrary curve. With the assumption 

of isentropic flow, so that i = Si, Eq. (6) can be differentiated to obtain 

p dq dq 


i - cu ^ + o/'r— - W^ 

p dq dq dq dq dq 


(7) 


Also the directional derivative of the pressure along this curve is 

dq dr dq d0 dq dz dq ' 

Evaluating the partial derivatives of the pressure in this equation by 


using Eq. (l), substituting this in Eq(7), and solving for S£, gives 


dq/ 


dW 

35 


2 , , TI (W a + cur) 2 

cu r + 1 dWr v 9 ' 


W 


W dt 


rW 


dr 1 ^ rW e + car2 ) de . 1 dWz dz 
dq W 


dt 


dq + W dt dq 


1 dh, 
+ — 


W dq 


£ <£ 3A /q\ 

" W dq ' ; 


Note that 


W r = sin a 


W z = Nm cos a 
W m = W cos (3 
Wg = W sin p 


and that 


da _ dq dm _ Win 
dt dm dt r c 

This can be used in evaluating the derivatives on the right side of Eq. (9) 
When this is done, we have 

I 

cu ^ (10) 


I 

§ = a * + b^+c^ + i^i 

dq dq dq dq W dq 


dq 
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where 


p p ,gtj 

= W cos a cos^p . w_sln_3 + sin a cos „ ■» . ^ sin p 

r„ r dm 


, W sin a cos 2 p ___ „ dW m 

h = - ^ + cos a cos B - — 

r c dm 

/ dW fl % 

c = W sin a sin p cos p + r cos p + 2co sin a 

V dm 


> (ID 


In equation (10) the values of the parameters h^ and A associated with a 
point inside the rotor is the value of that parameter at the inlet for the 
streamline which passes through the point. Then dh[/dq refers to the total 
enthalpy at the inlet as a function of the distance along the arbitrary curve 
q at the point considered. Further details of the derivation of Eq. (10) are 
given in [ 7] . 

For checking continuity of flow, it is necessary to calculate the density. 
Losses can be approximated by assuming a loss in relative total pressure Ap". 


With this 



where 


and 


_T_ _ W 2 + 2qjA - co 2 r 2 
Ti = ‘ 2c p T» 



2mA - pp 2 r 2 
2 c p t ! 


(14) 


Weight flow across a quasi -orthogonal can now be computed by 

'll 


w = 


H o 


pw n Andq 


(15) 
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where An is the distance between two stream surfaces on either side of the 
stream surface being analyzed, w is the weight flow between these surfaces 
per blade, and W n is the component of W normal to the surface defined by 
the quasi -orthogonal and the direction in which An is measured. The limits 
q _0 and q-j_ are the intersection of the quasi -orthogonal with the boundary of 
the flow passage. The flow area considered is indicated in Fig. 1. 

NUMERICAL TECHNIQUES AND PROCEDURE 

The procedure for utilizing Eq. (10) to obtain a solution on a prescribed 
stream surface is to estimate the parameters in Eq. (ll) so that Eq. (10) can 
be integrated with the condition that the continuity Eq. (15) be satisfied. At 
each quasi -orthogonal this determines a weight flow distribution, which can be 
used to obtain an improved streamline pattern. This leads to a better estimate 
of the parameters in Eq. (ll), and by iteration a flow distribution is obtained 
on the stream surface. 

Calculation of Parameters in Equation (10) 

The first step in the numerical integration of equation (10) is the 

dW dW fl 

numerical evaluation of the parameters a, 3, r„ , —Hi, and _ — for use in 

dm dm 

Eq . (ll). In order to evaluate the parameters a, 3, and r Q , a streamline 
geometry must be established. For this, fixed lines between flow boundaries are 
specified on the assumed stream surfaces, at several stations from inlet to 
outlet (for either a hub-to-shroud or a blade -to -blade stream surface). These 
Hines are the quasi -orthogonals along which Eq. (10) will be integrated. For 
an initial approximation to the streamlines, each quasi-orthogonal can be 
divided into a number of equal spaces. The quasi-orthogonais and initial stream- 
line assumptions are shown in Fig. 3 for a hub-to-shroud stream surface of a radial 
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inlet gas turbine. The success of the method is based on the fact that, for 
a reasonable assumed streamline pattern, the geometrical streamline parameters 
involved are not too different from those of the final solution. 

Along each assumed streamline we have r as a function of z at a number 

of points. It is desired to approximate — and £ from these values. For 

dz . dz 2 

any functions given at discrete points only, there are several ways a curve can 
be fitted through these values so as to approximate the original functions. A 
convenient method that has received much attention lately is the piecewise cubic, 
with matching first and second derivatives, usually referred to as a cubic spline 
function. Since for small slopes, the second derivative approximates the curvature 
of a function, the strain energy of a spline can be approximately minimized by 
minimizing ^ ^f "(x)j 2 dx, where f(x) denotes the curve described by the 
spline. The cubic spline has this property, as is proven in [9]. Thus the cubic 
spline is the mathematical expression of the shape taken by an idealized spline 
passing through the given points. In [9] a simple procedure is outlined for 
determining the spline fit curve when the coordinates of the points are given 
together with two arbitrary end conditions. One end condition that was used was 
that the second derivative at an end point is one half the second derivative at 
the next point. This is equivalent to bending the spline slightly beyond the 
last point, instead of just letting it be straight. An alternate end condition 
is the slope at an end point. The cubic spline function provided a simple 
analytical method of determining many of the parameters in the equations. The 
cubic spline was used to determine first and second derivatives, curvature, 
interpolated function values, interpolated derivatives, and for integration. 
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One further point concerning the spline fit should he mentioned; that is, 
the approximation to an actual spline curve is dependent on the slope not being 
too large. Experimentally, good results are obtained if the absolute value of 
the slope is not greater than one. In applying this method to streamlines on a 
radial turbine, there is a problem since the angle may be around -90° at the 
inlet. This is easily overcome by rotating the coordinate axes 45° so that the 
maximum slope is about one. 

After — and have been determined, a and — are given by 


dz 


dz 


a = tan 


-1 dr 
dz 

d 2 r 

dz 2 






(16) 


(17) 


dr 

dz 


The curvature rather than the radius of curvature is computed to avoid 

division by zero in case £ = 0. 

. 2 
dz 

The next quantity to be calculated is (3. Since 9 is known at each 
quasi -orthogonal along a streamline, the cubic spline can be used to approxi- 


mate 


d 0 

The angle (3 is then calculated from 

3 = tan' 




(18) 


For an initial calculation, W may be assumed constant throughout the 

dW„ dW 0 

rotor. Since Wm = ¥ cos |3 and Wa = W sin (3, the quantities - — and - — 

a dm dm 

can now be approximated using the cubic spline. Now all quantities necessary 
for the integration of Eq. (10) are determined. 
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Integration of Equation (10) 

Given an initial value of W, Eq. (10) can be accurately integrated by 
means of a Runge-Kutta method. Reference [7] gives further details on the 
method. Since the initial value is not known, a guess of some kind must be 
made. This results in the velocity distribution along the quasi -orthogonal. 
Equations (12) to (14) can be used to compute the integrand in Eq. (l5), which 
can then be numerically integrated using the cubic spline approximation. The 
computed total weight flow is then compared with the actual weight flow. If 
the computed weight flow is too small, the initial velocity is increased, and 
vice versa. Then the velocity distribution and the weight flow are recalculated. 

A few iterations will determine the initial velocity that will give the correct 
weight flow. 

Calculation of Improved Streamlines 

From Eq. (15) the weight flow distribution along the quasi -orthogonal can 
also be obtained. Along any quasi -orthogonal, inverse interpolation (by a spline 
approximation) can be used to determine the streamline spacing that will give 
equal weight flow between any two adjacent streamlines (see Fig. 4). When this 
is done for every quasi -orthogonal from inlet to outlet, a new estimate for the 
streamline pattern is obtained. However, if this streamline pattern is used for 
further iterations the procedure does not converge due to the fact that the 
streamlines become rough and uneven. This may be avoided by using only a fraction 
of the calculated streamline corrections, perhaps one-tenth or less. With this, 
the streamlines remain smooth, and a solution is reached in a single computer 


run. 
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A Quasi -Three -Dimensional Solution 

The procedure described gives a solution for any known stream surface. 

Since the stream surfaces are not known they must be approximated. For this 

reason a true three-dimensional solution is not obtained, but only what is 

usually referred to as a quasi-three-dimensional solution. One procedure 

that can be used is as follows. Assume that the mean hub-to-shroud stream 

surface has the same shape as the mean blade surface. A deviation from this 

may be specified near the inlet if there is an incidence angle. For this 

initial approximation. An in Eq. (15) is equal to - t fl . With these 

n 

assumptions the QO Method is used to obtain a meridional solution. The 
streamlines obtained from this solution are used to define surfaces of 
revolution from blade -to -blade, with An (in Eq. (15)) for these blade-to- 
blade surfaces determined by the streamline spacing in the meridional solution. 

A separate blade -to -blade solution could be obtained for each streamline; 
however, in most cases three blade -to -blade solutions would be adequate, one 
at the hub, one at the mean streamline between hub and shroud, and one at the 
shroud. This gives a relative velocity distribution throughout the passage. 
Further details of the procedure outlined above are given in [7 and 8j. 

The solution could be refined by using the blade -to-blade solution to 
obtain a better approximation to several hub-to-shroud stream surfaces, con- 
tinuing until the hub-to-shroud and blade -to-blade solutions were completely 
consistent. This has not been done. 

NUMERICAL EXAMPLE 

The method outlined has been applied to the analysis of a small radial 
flow turbine. The calculations were made on a digital computer. The hub-shroud 
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profile and quasi -orthogonals in the meridional plane are shown in Fig. 3, 
together with the equally spaced streamlines used for the initial assumption. 
The blade has radial elements, except near the trailing edge of the splitter 
blade, where the taper is not the same on both sides of the blade. The 
pertinent data for the case analyzed is given in the following table. 


DATA FOR NUMERICAL EXAMPLE 


Total number of blades 

22 

Number of splitter blades 

11 

Tip radius, r^, in. 

3.01 

Rotational speed, r.p.m. 

38 500 

Fluid 

Argon 

Weight flow, lb/sec 

.611 

Inlet total temperature Tj_, °R 

1950 

Absolute tangential velocity at inlet, Vq j_, ft/sec 

883 

Inlet relative flow angle, Pj_, degrees 

-25.1 

Inlet total pressure, pj, psia 

12.97 

Loss of relative total pressure between inlet 

.25 

and outlet, Ap", psi (Assumed to vary linearly 


along streamline . ) 



The normal blade thickness was given by means of tabulated values on a grid. 
Blade thickness at any given point was obtained by linear interpolation. It 
was assumed that h{ and A were constant across the inlet. 

At the inlet, the hub-shroud stream surface was assumed to deviate from 
the blade surface in order to agree with the flow direction coming into the 
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rotor. Let r^ denote the radius where the mean stream surface is assumed to 
deviate from the mean "blade shape. Equation (13) of [5] gives an approximate 
equation for determining r^, which may be written as follows: 


r b = r i e 


-0.71 A 0 


(19) 


The stream surface was assumed to vary cubically with r between r^ and r^. 

The first step in the solution is to obtain a meridional plane solution. 

This solution was extended downstream of the blades, by assuming a stream 
surface extending downstream from the blades at an angle determined by the 
blade angle corrected for blockage. Twenty meridional streamlines were calculated 
by the Q0 Method. The meridional streamline pattern obtained at hub, mean, and 
shroud is shown in Fig. 5. This gives the coordinates for a blade -to -blade 
surface at the mean meridional streamline, as well as the streamline spacing 
on the quasi-orthogonals at hub, mean, and shroud. With this, the normal 
streamline spacing and the blade coordinates were calculated at each of the three 
surfaces, and three blade -to-blade solutions were obtained by the Q0 Method. 

The three blade -to-blade solutions were extended upstream of the inlet, because 
of nonuniform conditions at the inlet. The solution gives the Q and m 
coordinates of the blade -to-blade streamlines, which are plotted in Fig. 6. 

The complete quasi-three-dimensional solution was obtained in a single computer 
run, taking 6 minutes on the IBM 7094 computer. 

Figure 7 shows the blade loading on main and splitter blades at hub, mean, 
and shroud. The velocities are increasing over most of the blade. However, 


there is a large negative velocity gradient on the pressure surface at the inlet 
near the hub. It is even more severe at the end of the splitter blade near the 
shroud. This indicates the splitter should be extended further from the aero- 
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dynamics point of view. However, the splitter is already longer than needed 
at the huh, and stress consideration may not permit making the splitter extend 
beyond the point where it ends on the hub. Figure 8 shows the absolute velocity 
distribution on the stationary shroud. Here, there are large negative velocity 
gradients. Of course, a negative velocity gradient cannot be eliminated 
entirely with the leaving absolute velocity lower than the inlet velocity, 
although it can be minimized by careful design. On the pressure surface the 
nearly linear distribution is about the best possible for avoiding flow separa- 
tion, under these circumstances. However, along the shroud there is an increase, 
and then a greater decrease in a shorter distance. It is difficult to avoid 
something like this if there is to be any loading of the blades at the shroud. 

This type of analysis could be very useful as a design tool, since it 
points up the location of possible flow separation. Modifications can be made 
in the geometry to improve on the velocity distributions until a good design 
is evolved. 

CONCLUDING REMARKS 

A method of analysis of turbomachines is summarized that is suitable for 
computer programing. The method, which has been reported in references 7 and 8, 
is based on a streamline analyisis of a stream surface using quasi-orthogonals. 

A quasi-three-dimensional solution can be obtained by using the method first for 
a hub-to-shroud analysis, followed by blade -to-blade analyses at hub, mean, and 
shroud. As an example, the method was applied to a radial inflow turbine with 
splitter blades. The complete quasi-three -dimensional solution was obtained in 
a single computer run. The FORTRAN computer program used for the numerical 
example may be obtained from the author by anyone who is interested in using the 


method. 
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Figure 7. - Blade velocity distribution for numerical example. 







Figure 8. - Shroud absolute velocity distribution for numerical 
example. 
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